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Abstract 

We consider the scaling limit of a generic ferromagnetic system with a continuous phase 
transition, on the half plane with boundary conditions leading to the equilibrium of two 
different phases below criticality. We use general properties of low energy two-dimensional 
field theory to determine exact asymptotics of the magnetization profile pcrpcrdicularly to 
the boundary, to show the presence of an interface with endpoints pinned to the boundary, 
and to determine its passage probability. The midpoint average distance of the interface 
from the boundary grows as the square root of the distance between the endpoints, unless 
the reflection amplitude of the bulk excitations on the boundary possesses a stable bound 
state pole. The contact angle of the phenomenological wetting theory is exactly related to 
the location of this pole. Results available from the lattice solution of the Ising model are 
recovered as a particular case. 



1 Introduction 



Interfacial phenomena at boundaries are a subject of relevant interest for both theory and 
applications. On the theoretical side, the one this paper is concerned with, the effects of the 
boundary on an interface separating different phases of a statistical system have been extensively 
studied using phenomenological, mean field, renormalization group and other approximation 
methods ([1-8] is a certainly incomplete list of review articles). The only exact result that has 
been available concerns the Ising model on the half plane [HI [10], a circumstance that, while 
confirming a specificity of the two-dimensional case, raises the question about the role of Ising 
solvability in these exact findings. 

We show in this paper that exact results including those of [10] as a particular case are 
obtained quite generally for any two-dimensional model exhibiting a continuous phase transition. 
This is done extending to the half plane the non-perturbative field theoretical approach recently 
used in [11] to study phase separation on the whole plane. As in that case, general exact results 
emerge because, when its end-to-end distance R is much larger than the correlation length, the 
interface is described by a single particle (domain wall) state, in a low energy limit leading to a 
general solution. In this way, the fluctuations of the interface turn out to be ruled by the low 
energy singularity of the matrix element of the order parameter field (as for the whole plane), 
with the fields pinning the interface endpoints to the boundary producing boundary reflection 
and an average midpoint distance from the boundary of order y/R. 

The result changes qualitatively if boundary and domain wall excitation admit a stable bound 
state, which becomes dominant in the spectral sum at low energies and bounds the interface to 
the boundary. The contact angle and the spreading coefficient of the phenomenological theory 
of wetting then emerge in a completely natural way within the field theoretical formalism. 

The paper is organized as follows. In the next section we illustrate the field theoretical 
setting and derive the results for the unbound interface. Section 3 is then devoted to the effects 
produced by the bound state and to the characterization of the wetting transition, while section 4 
contains some final remarks. 

2 Interfaces on the half plane 

Consider a ferromagnetic spin model of two-dimensional classical statistical mechanics in which 
spins take discrete values labelled by an index a = 1,2, ... ,n. The energy of the system is 
invariant under global transformations of the spins according to a symmetry whose spontaneous 
breaking below a critical temperature T c is responsible for the presence on the infinite plane of 
n translation invariant pure phases; we denote (• • • ) a statistical averages in the phase a. 

Assuming a continuous transition, we consider the scaling limit below T c , corresponding to 
a Euclidean field theory defined on the plane with coordinates (x,y), which can be seen as the 
analytic continuation to imaginary time of a (l+l)-dimensional relativistic field theory with 
space coordinate x and time coordinate t = iy. If H and P are the Hamiltonian and momentum 
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Figure 1: Elastic scattering (reflection) of a kink off the boundary (a), and interface pinned at 
the boundary (b). 



operators and $ a field of the theory, translation invariance on the plane yields the relation 

$(x, y) = e ixP+yH ^(0, 0)e~ ixP ~ yH . (1) 

The (l+l)-dimensional theory possesses degenerate vacua |0) a associated to the pure phases of 
the system. The elementary excitations correspond to stable kink states \K a b(9)) interpolating 
between different vacua |0) a and |0}&. We introduced the rapidity variable 9 which conveniently 
parameterizes the energy and momentum of the kinks as (E,p) = (mcosh#,msinh#), m being 
the kink mass or inverse correlation length. The trajectory of the kink on the Euclidean plane 
corresponds to a domain wall between the phases a and b. Multi-kink excitations take the form 
\K aai (6i)K aia2 (92) ■ ■ ■ K an _ 1 b(9 n )) . Within the scattering framework [12] we consider, these are 
asymptotic states, incoming if considered long before the collisions among the kinks, outgoing 
if considered long after, and their energy is simply ^^ 1 mcosh^j. 

Consider now the system on the half-plane x > 0. We denote by B a a boundary condition 
at x = which is y-independent and breaks the symmetry of the bulk in the direction a in order 
parameter space; this can be realized applying a constant boundary magnetic field pointing in 
the direction a. We denote (• • • )s statistical averages in presence of the boundary condition 
B a . Preservation of translation invariance in the y direction yields energy conservation in the 
(1 + l)-dimensional picture. The bulk excitations are still the kink states described for the full 
plane case, but now they are restricted to x > 0; we indicate this restriction by a subscript 
B a . Hence |0)b o denotes the vacuum (no excitations in the bulk) on the half-plane with the 
boundary condition B a . If a is the spin field, the magnetization {cr(x,y))B a = B a (0|o"(x, y)\0)B a 
points in the direction a and depends only on the distance x from the boundary; in particular 

lim (a(x,y)) Ba = (a} a , (2) 

x— >oo 

where (a) a is the constant magnetization in phase a on the full plane. The state |0)s a is 
eigenstate of the Hamiltonian Hs a of the system on the half line. We consider the case in which 
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boundary conditions B a and Bb are related by the symmetry, so that |0}_b o and |0)b 6 have the 
same energy Eb- 

The asymptotic scattering state \Kb a {9)) B a corresponds to an incoming kink (travelling to- 
wards the boundary) if its momentum is negative, i.e. if 9 < 0. If its energy is lower than the 
energy 2m needed to produce two kinks upon interaction with the boundary, it will simply be 
reflected into an outgoing kinlfl with rapidity —9 (Fig. la). The state \K~b a {9))B a is eigenstate 
of HB a with eigenvalue Eb + mcosh#. 

We are now ready to set up the configuration we want to study, namely a boundary condition 
which is of type B a if |y| > R/2 and of type Bb if \y\ < R/2. The interest of such a boundary 
condition, that we denote B a j, a , is easily understood observing that the limit for x — > oo of the 
magnetization profile (cr(x, 0))B aba has to tend to (a) a if R is finite, and to (a)b if R is infinite. 
The natural way to account for this situation is to expect the formation of an interface pinned 
at R/2 and —R/2 on the boundary, separating an inner phase b from an outer phase a (Fig. lb), 
and whose average distance from the boundary at y = diverges with R. The remainder of this 
section is devoted to see how such a picture indeed emerges within our general field theoretical 
framework. 

Technically the change from the boundary condition B a to Bb at a point y is realized starting 
with B a and inserting on the boundary a field fi a b(0, y) which acting on the vacuum |0) s a creates 
kink states interpolating between phase a and phase b. Hence the simplest non-vanishing matrix 
element of the boundary field fi a b is 

B a (0M0,y)\K ba (9)) Ba =e-y mcoshe B a (0\f, a b(0,0)\Kba(0))B a ^ e^™ he ^(0) . (3) 
The partition function of the system with boundary condition B a b a reads 

Z = Ba (0\v ab (0,R/2)» ba (0,-R/W)B a = J £\^(9)\ 2 e- mRcoshd + 0(e- 2mR ), (4) 

where the last expression is obtained expanding over an intermediate set of outgoing kink states 
and retaining only the lightest (single kink) contribution which is leading in the large mR limit 
we will consider from now on. Since the above integral is dominated by small rapidities and J-^ 
is expected to behave as0 

F^9) = a9 + 0{9 2 ), (5) 

the partition function becomes 

/■oo jq I |2 —mR 

Z ~ \ a \ 2 / — 9 2 e —^( 1 + e2 /2) = _N£ (g) 

' ' Jo 2tt 2V2^ (mR) 3 / 2 V ' 



x As emphasized in [T3], the analogies between bulk and boundary scattering become evident thinking of the 
boundary as the propagation of an infinitely heavy particle sitting at x — 0. 

2 Linear behavior of matrix elements at small rapidities in two-dimensional theories is well known. Within the 
framework of integrable boundary field theory [13) exact examples can be found in [14]. More generally, see [15] 
about matrix elements in integrable theories. 
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The magnetization profile along the x axis is given by 

(a(x,0)) Baba = ± Ba (0\ t , ab (0,R/2)a(x,0) f i ba (0,-R/2)\0) Ba (7) 

i /"+00 jn in 

where in the last line we have taken mR ^> 1 to project on the one-kink intermediate states, but 
also mx 3> 1 to be able to treat cr(x, 0) as a bulk field which satisfies ([T|) and is evaluated on bulk 
kink states (whose rapidities take both positive and negative values). In other words, for mx 
large the only effect of the boundary on the magnetization comes from the boundary changing 
fields at (0, dzR/2); in their absence one would simply observe the constant value {cr) a . The bulk 
matrix element of the spin field between one-kink states is related by the crossing relational 

<tfab(0i)k(O, 0)\K ba (9 2 )) = + m - 9 2 ) + 2*5(01 - 2 ){a) a , (8) 

to the form factor 

" 2 ) = a ^\<y^MKab(ei)K ba {02)) ■ (9) 

As already observed in [11] for the case of phase separation on the whole plane, it is crucial that 
quite generally, due to non-locality of the kinks with respect to the spin field, F a {9) possesses 
an annihilation pole at 9 = in with residue [16] 



- i Res e=lw F a (e) = (a) a - {a) b = A (a) . (10) 
Since mR is large (|7|) is dominated by small rapidities and (0), and (flU]) lead to 

(a(x,0))B aba ~2(o) a + iA(a) ^e'™* £j ^ _M|_ e *-^-(^l)«] . (11 ) 

Differentiation removes the singularity of the integrand and gives 

\ a \2 e -rnR 

d mx (a(x,0)) Baba ~ -A(g) / 2 x 2z 9(x)g(-x) 



Ay/2 2 _ z 2 [2m~ 



y(*)^==z 2 e- z , Z = V~W X ( 12 ) 

VirmR V it 



where we used (© and 



Q2 //i I „-™~fl 



5 (x) = / d0 e -mM' = z e - 2 /2 . (13) 

Integrating (fl2j) with the asymptotic condition (a(oo,0)) Baba = (a) a gives 

0)) Baba ~ (a) 6 - A (a) f * e"* 2 - jT e""^ , mx»l. (14) 

3 Crossing a particle from the initial to the final state (or vice versa) involves reversing the sign of its energy 
and momentum [12], namely an in rapidity shift. The delta function term in |[8} is a disconnected part arising 
from annihilation of the two kinks. 
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From this result we can compute exactly limR_ s . 00 (a(^ (mR) s ,0)) B aba , obtaining (a)b for < 5 < 
1/2, (a) a for 5 > 1/2, and the r.h.s. of flU} with z = ay/2 for 5=1/2. For (a) a = -(a) b = {a) + 
these are precisely the limits obtained from the lattice in [17\ [9] for the Ising model on the half 
plane with boundary spins fixed to be positive for \y\ > R/2 and negative for \y\ < R/2. 

The derivative (|12p of the magnetization profile is peaked around z = 1, confirming the 
presence of an interface whose average distance from the boundary increases as \J R/m. It is also 
easy to see that the result for the magnetization profile is consistent with a simple probabilistic 
interpretation. Since we are computing the magnetization on a scale R much larger than the 
correlation length and far away from the boundary, we can think of the interface as a sharp 
separation between pure phaseqj, and write 

rx poo 

(a(x,0))B aba ~ {c) a I dup(u) + (a) b dup(u), mx > 1 , (15) 

JO Jx 

where p(u)du is the probability that the interface intersects the rr-axis in the interval (u, u + 
du), so that the two integrals are the left and right passage probabilities with respect to x. 
Differentiating and comparing with (|12p gives the passage probability density 



V^ 2 e- 2 , (16) 



which correctly satisfies / °° dx p{x) = 1. 



3 Wetting transition 

The results of the previous section are modified if the kink-boundary system associated to the 
asymptotic state \K a b(9)) s b admits a stable bound state |0)b^, corresponding to the binding 
of the kink K a b on the boundary Bb- As usual for stable bound states [12], such a binding 
will correspond to a "virtual" value 6*o of the kink rapidity leading to a bound state energy 
Eb + mcosh#o real and smaller than the unbinding energy Eb + fn. This amounts to taking 
9q = iu with < u < 7r, so that 

Eb> = Eb + to cos u . (17) 

The existence of the bound state manifests in particular through a simple pole in the elastic 
scattering amplitude of the kink off the boundary, which reads TZ(6) ~ ig 2 /(0 — iu) for 6 — > iu, 
with g a kink-boundary coupling constant (Fig. 2a). This pole is inherited by the matrix element 
([3]), for which we hav|] (Fig. 2b) 

M e ) = B a (0\tiab(0,0)\K ba (0)) Ba ~ 7 ^—B a (0\Vab(0,0)\0} B > , 9 ^ iu . (18) 



4 It has been shown in [TT| how the internal structure of the interface arises from subleading terms in the large 
mR expansion. 

Exact solutions exhibiting boundary bound states poles can be found in [13] for scattering amplitudes and in 
[14] for matrix elements. 
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Figure 2: The boundary bound state (double line) originating in kink-boundary scattering (a), 
and a pictorial representation of equation ([15]) (b). 

The boundary bound state affects the results of the previous section for the boundary condi- 
tion B a b a because the leading low-energy contribution in the expansion over intermediate states 
now comes from |0}b' rather than from \K\, a {9)) Ba . So the partition function becomes 



Z = Ba (0\ fJ , ab (0,R/2)» ba (0,-R/2)\0) Ba = \ Ba (0M0,0)\0) B ,\ 2 e- mRcosu + O(e- mR ), (19) 



We see then that, as a consequence of ([2]), the magnetization profile now tends to {cr) a at large 
mx, in contrast to what obtained in the previous section, where it tended to {a)b for R large 
enough. This corresponds to the fact that now the asymptotic behavior is determined by the 
state in which the interface, and then the phase b, are bound to the boundary, while before the 
dominant state was that in which phase b extended to an average midpoint distance of order 
\f~R from the boundary. 

Consistency of the asymptotic expansion requires that the corrections to (|20p vanish as 
R — > oo. For mx large, the first of these corrections is that due to the \K\j a {0)) intermediate 
states given in (J7]). The Z in the denominator, however, is now (I19D rather than ([6]), so that 
the correction behaves as e mF! -( cosu - 1 ) a t large R. Hence, if u approaches 0, i.e. if the interface 
approaches the unbinding point, consistency requires that R diverges faster than 1/u 2 . If we 
adopt a vocabulary within which b is a liquid phase and a a vapor phase, we can say that as 
u — > a thin layer of the liquid phase spreads all over the boundary. 

The relation with the usual characterization of interfacial phenomena at boundaries be- 
comes more transparent if we consider the situation usually referred to as "partial wetting", 
corresponding to a drop of liquid sorrounded by a thin layer of liquid adsorbed on the rest of 



and the magnetization profile 




(20) 
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Figure 3: Splitting and recombination of the boundary bound state B' a corresponds to "partial 
wetting" , in which a drop of phase b makes an equilibrium contact angle 9 e with the boundary. 
Equation (|17p with u = 9 e gives the surface tension balance condition at the contact points. 

the boundary (see e.g. [8]). In our formalism this amounts to splitting and recombination of the 
boundary bound state B' a (Fig. 3). Considering that the kink mass m is the surface tension of 
the interface [TT], that Eb is the surface tension between the boundary and the drop, and that 
Egi is the surface tension between the boundary away from the drop and phase a, we recognize 
in (|17p the Young equilibrium condition at contact points (see e.g. [1] and references therein), 
with u playing the role of the equilibrium contact angle 6 e (Fig. 3). In addition, the combination 
m(cosu — 1) encountered a moment ago is recognized as the so called "equilibrium spreading 
coefficient" (see [8]). We also see that interface unbinding at u = corresponds to vanishing of 
the contact angle, namely to the usual characterization of the wetting transition point (passage 
from partial to complete wetting). 

The boundary bound state is a property of the theory with translationally invariant boundary 
condition B^. Parameters of this theory are the temperature, related to the kink mass as 
m oc (T c — T) u , and a coupling A entering the boundary term A J dy(j)(0,y) of the classical 
reduced Hamiltonian. If X is the scaling dimensiorJ§ of the boundary field <f>(0,y), u is function 
of the dimensionless combination \/m}~ x . If A is kept fixed, the condition u = determines a 
wetting transition temperature T w (\) < T c . 

The results ©, (|19p and (|20p account for those reported in [9l[T0] for the particular case of 

an Ising model with boundary condition B^ |_ and coupling between the boundary spins and 

their nearest neighbors different from the coupling within the rest of the lattice; this modified 
coupling corresponds to the boundary parameter A in this case. The generality of our results 
also explains why approximated treatments of other models resulted in findings similar to the 
Ising ones (see |10] and references therein). 

6 The exponents v and X are known exactly from bulk [18] and boundary [19] conformal field theory, respec- 
tively. 
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4 Conclusion 



In this paper we studied the scaling limit of a generic ferromagnetic system with a continuous 
phase transition, below criticality and on the half plane, with boundary conditions favoring one of 
the phases along an interval of length R, and a different phase outside this interval. We used field 
theory to determine exact large R asymptotics of the magnetization profile perperdicularly to 
the boundary at the middle of the interval. We showed that, generically, the large R asymptotic 
behavior corresponds to the presence of an interface pinned at the boundary condition changing 
points, with an average midpoint distance from the boundary which grows as y/R. The passage 
probability density of the interface has the gaussian form found in for the whole plane, 
modified by a quadratic factor which accounts for the presence of the boundary. These results 
are modified if the scattering on the boundary admits a stable bound state, which then becomes 
leading at low energies and corresponds to the binding of the interface to the boundary. In 
this case we showed how field theory accounts at a fundamental level for the contact angle and 
spreading coefficient of the phenomenological wetting theory. 

These results follow from general low energy properties of two-dimensional field theory. In 
particular, the annihilation singularity of the spin field matrix element on one-kink states and 
the boundary-kink bound state pole play a key role in determining the asymptotics of the 
magnetization profile in the unbound and bound regimes, respectively. 

Additional interfacial properties, such as the internal structure arising from subleading terms 
of the large R expansion or double interfaces appearing in some models for particular choices of 
boundary conditions, can be analyzed in the same way it was done in [11] on the whole plane; 
we refer the reader to that paper on these points. 
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